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s ! ■ Abstract. In these lectures I introduce the basics of HERA physics and give a survey 

(— i ' of the major aspects, discussing in somewhat more depth the subject of low x physics. 

1—5 ■ PACS: not given 

m '■ 

j> \ 1 Introduction 

7— I '. 

The study of deep inelastic lepton-proton scattering has produced some of the 
major underpinnings of the Standard Model. For example, the quark-parton 
model took shape in the light of the deep inelastic scattering (DIS) experi- 
£S) • ments jlj begun at SLAC in the late 1960s. Going even further back, the scat- 

' tering of energetic "simple" a particles from the nuclei in a thin gold foil, carried 

out by Geiger and Marsden in Manchester in 1909, led to the concept of the nu- 

| clear atom || and is clearly analogous to deep inelastic lepton scattering in 

modern particle physics. 

' In my lectures I first gave an overview before covering one particular area 

in more detail. In the interests of producing a more coherent write-up, I have 
changed the original order. I first outline the HERA accelerator and detectors 
before discussing the theoretical techniques used to derive the basic formulae 
used in the study of electron-quark scattering. Then I give an overview of the 
main areas of HERA physics, going into somewhat more detail in one particular 
area, that of electron-proton scattering when the interacting quark has a low 
fraction of the original protons momentum, so-called "low- a;" physics, and re- 
lated areas such as diffraction. Here, recent theoretical developments and model 
building are changing our perception and improving our understanding of the 
very rich phenomenology arising from the many different but related channels 
that can be explored at HERA. I conclude with a discussion of the recent up- 
grades to both the accelerator and experiments in the HERA II programme and 
the main areas of physics that they will address. 



X 



I Delivered at 7th Hellenic School on Elementary Particle Physics, Corfu Summer Institute, 
September 2001 
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2 Introduction to the HERA machine and experiments 

HERA is a unique facility, colliding beams of electrons or positrons with pro- 
tons at high energy. The protons are accelerated and stored in a ring of super- 
conducting magnets; until 1998 the protons were accelerated to 820 GeV and 
subsequently to 920 GeV. The electron or positron ring is normal conducting 
and beams are stored at 27.5 GeV. HERA began operation in 1992 and con- 
tinuously improved its performance in successive years, as illustrated in Fig. [j], 
which shows the luminosity delivered to the ZEUS experiment in each year of 
running. Because of difficulties with the electron lifetime, the great majority of 
data has been taken with positron beams; only 32 pb _1 of electron data has 
been collected by ZEUS compared to a total of 148 pb _1 of positron data. 
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Fig. 1. The luminosity delivered to the ZEUS detector as a function of days of 
the run, shown separately for each year of running. A continuous improvement 
in the performance of HERA is evident, except for 1998, in which electrons were 
used rather than positrons. 

There are two "general-purpose" detectors at HERA, HI and ZEUS. In ad- 
dition the HERMES experiment uses a gas target to examine polarised electron 
or positron-polarised proton scattering, and the HERA-B detector, designed to 
study CP violation in the B sector. The latter two are not discussed further in 
these lectures due to lack of time. Both HI and ZEUS have a rather similar con- 
figuration, as far as possible enclosing the full solid angle with tracking detectors 
surrounded by calorimetry. Because of the large asymmetry between the proton 
and positron beam energies, the energy flow is predominantly in the proton, or 
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"forward" direction, so that the detectors are asymmetric, with thicker calorime- 
tery and a higher density of tracking detectors in the forward direction. Figure ^| 
shows a diagram of the ZEUS detector, illustrating that its general structure, 
with the exception of the more complex forward instrumentation, is very typical 
of modern Air detectors, such as those at LEP and the Tevatron. 




Fig. 2. A vertical section through the ZEUS detector. The black cylinders in 
the very centre represent the silicon microvertex detector, which was installed in 
2000. Surrounding this are the tracking detectors (CTD, FDET and RTD) and 
surrounding these are the uranium/scintillator calorimeters (FCAL, RCAL and 
BCAL). These are surrounded by an iron-scintillator backing calorimeter (BAC) 
which also acts as a flux return for the superconducting solenoid which surrounds 
the CTD. Muon chambers (FMUON, with associated toroidal magnets, BMUO 
and RMUO), a Veto Wall to veto off-momentum protons and a concrete shield 
complete the detector. 

With the advent of HERA, the accessible phase space in the kinematic in- 
variants Q 2 (the virtuality of the exchanged virtual photon) and x increased 
by approximately three orders of magnitude in each variable compared to what 
was available at earlier fixed-target experiments (see Fig. ||). This extension in 
kinematic range has opened up qualitatively new fields of study, both at high 
and low Q 2 . 

Since HERA is a high-energy lepton-hadron collider, it clearly gives access 
to hard processes in both the strong and electroweak interactions. HERA is very 
sensitive to the production of new particles that can be formed by the fusion 
of leptons and quarks, e.g. leptoquarks or many of the particles predicted by 
R-parity-violating supersymmetry. It is also very sensitive to any small changes 

http : / /link . springer .de/link/ service/ journals/10105/ index.html 




Fig. 3. The kinematic plane in x and Q 2 for experiments probing the parton 
distribution of the proton. The regions explored by each experiment are shown 
in a variety of shadings as shown in the legend. Hadron-hadron collisions are also 
able to measure the proton structure, predominantly at high x and high Q 2 . 



in the pattern of the electroweak interaction predicted by the Standard Model. 
Since the colliding leptons are point-like, HERA allows complete control over 
the conditions of the collision by varying the Q 2 , and thereby the size, of the 
probe. At high Q 2 , HERA is a probe of the complex structure of the proton 
via the point-like coupling of the photon. At low Q 2 , the photon becomes large 
and evolves its own complex structure which can be probed using the point- 
like interactions between its parton constituents and those of the proton. The 
hadronic nature of the photon under these conditions gives rise to hadron-hadron 
interactions with large cross sections; by also analysing diffractive interactions, 
in which the proton can be violently struck but remain intact, the rich structure 
and phenomenology of the strong interactions can be explored. 

One thing that HERA physics cannot do is be simple. Unlike the situation 
in electron-positron annihilation, energy and quantum numbers are transferred 
between the colliding particles, each of which has its own conserved quantum 
numbers of lepton and baryon number. This means that the single annihilation 
energy necessary to describe most of electron-positron or high-energy hadron- 
hadron collisions is insufficient; two invariants are required. These can be picked 
from several different possibilities, the most common of which are x and Q 2 . 
Others include W 2 , the square of the energy of the hadronic final state, s, the 
squared centre-of-mass energy of the electron-proton system, or y, the inelas- 
ticity, which in the rest frame of the proton is the energy transferred from the 
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electron to the proton. Only two of these variables are independent; their defini- 
tions are given in Eq. 1 in terms of the initial- and final-state four-vectors of the 
electron, k and k' respectively, and the same quantities, p and p', for the proton. 
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s = (p + k) 2 

f = - q 2 = -(k'~k) 2 



p ■ k 

w 2 = ( P ') 2 = {p + q? 

Energy-momentum conservation implies that 

Q 2 

x = ^ — 

2p-q 

so that, ignoring the masses of the lepton and proton: 



(2) 



t1 



W = Q 2 (4) 

x 

Another thing that one cannot expect from HERA is to discover the Higgs. 
HERA physics is dominated by the collisions of the lightest quarks with the 
lightest leptons, so that it is the worst place to look for new particles whose 
couplings are proportional to mass. Neither can HERA produce new particles 
with mass close to the centre-of-mass energy unless they have the quantum 
numbers of leptoquarks; a large proportion of the energy must be carried off by 
the final-state lepton and baryon in the t channel and is not available for the 
production of new particles. 

Finally, it cannot be expected that HERA physics will be simple. Many of the 
simplest problems in QCD have already been studied in detail. The remaining 
ones, such as the problem of confinement, are of the highest importance but 
extremely difficult to study in practice, since they manifest themselves in regimes 
in which the strong interaction is really strong. 



3 The formalism of Deep Inelastic Scattering 

It is instructive to remind oneself that the basic formalism of DIS can be rel- 
atively easily derived^ from the QED treatment of spin-i - spin-i scattering. 
It is convenient to work in a frame in which the interaction of the virtual pho- 
ton with the constituents in the proton can be considered as incoherent, i.e. 
the characteristic time of the "/*q interactions is much shorter than any interac- 
tions between the partons. A suitable frame is the infinite-momentum frame of 
the proton, which at HERA can be approximated by the lab. or centre-of-mass 
frame. In such a frame, Lorentz contraction reduces the proton to a "pancake" 

1 In my lecturers I spent considerable time in deriving many of the standard DIS equations 
from first principles. My approach was based on the use of the Mandelstam variables and 
follows closely the treatment in Chapters 6 - 10 of Halzen and Martin M. Only an outline of 
the main points is reproduced here. 
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and time dilation increases the lifetimes of the fluctuating partons so that the 
proton constituents are effectively "frozen" . Provided that the quarks have neg- 
ligible effective mass, i.e. have small rest mass and are asymptotically free, and 
that Q 2 <C fcy, then the interactions can be considered incoherent ||. 

The basic process of an electron scattering incoherently from a quark is now 
identical to the classic QED calculation of electron-muon scattering. Writing the 
initial and final four-momenta of the electron as k,k' and those of the proton as 
p,p', the standard Feyman rules allow us to write the matrix element as 

M. = —e 2 u(k')^u(k)—u(p')^ u u{p) (5) 
1 

Using standard trace techniques and ignoring mass terms leads to 

\M? = j^rjy [(k'-p')(k-p) + [(k'-p)(k-p')\ (6) 

It is very convenient to use the Mandelstam variables, s, t and u, since not only 
can they easily be evaluated in any frame, but also there are several useful 
relations between them and the more usual DIS variables of Eq. 1 that simplify 
the algebra. The Mandelstam variables are defined in terms of the four-vectors 
as: 

s = (k+p) 2 ~ 2fc-p~ 2k' -p', (7) 
* = {k~k'f ~ -2k -k' ~ -2p-p', (8) 
u = (k — p 

so that Eq. || simplifies to become 

\M\ 2 ^{s 2 + u 2 ) (10) 

Note that crossing, the replacement of s by t and vice-versa, leads to the well- 
known formula for e + e~ — > /i + . 

The matrix element can be converted to a cross section by using the standard 
formula for 2^2 scattering, 

to give 



u = (k - p'f ~ -2k ■ p' ~ -2k' ■ p, (9) 



da • 1 



' (§ 2 +u Z ), (12) 



dt 8irs 2 t 2 

where the hatted variables represent the Mandelstam variables for the subprocess 
in question, which for e — /x scattering are identical to the unhatted variables, but 
which for eq scattering are not. Considering now eq scattering, the Mandelstam 
variables satisfy 

s + £+m = 0, (13) 



http : //link. springer . de/link/service/ journals/ 10105/ index. html 



EPJdirect Al, 1-11 (2008) Springer- Verlag 



8 



which, assuming that the quark brings a fraction x of the protons energy into 
the eq collision can be written as 

x(s + u)+t = 0. (14) 

Equation [l^ can be converted into a double-differential cross section using the 
appropriate S function to give 

which, using the appropriate Jacobian, can be written in terms of s,t and u as 
d 2 a 2~Ka 2 xe 2 

— = -^(s> + u>W + X (s + u)), (16) 

where a is the fine-structure constant, a — e 2 /Aw and e q is the charge of the 
struck quark in units of the electron charge. We can then write the total ep cross 
section as the incoherent sum of all possible eq scatters, i.e. 

/ cPtr \ f . , , f cPa \ 



\dtdu 



e.X 



where fi(x) is the density distribution inside the proton of quark i between x 
and x + dx. 

We can also treat the overall inclusive DIS process ep — ► eX from first prin- 
ciples using the Feynman rules provided that we take cognizance of the fact that 
the proton is not a point-like particle. To do that, we parameterise the proton 
vertex contribution to the matrix element in the most general way possible in 
terms of a hadronic tensor, W^ v , given by 

W, a ,= (-g ltv + Fl {x,Q 2 ) + P -^F 2 {x,Q 2 )~ ie^ pt F 3 (x,Q 2 ) 

(18) 

where 

Pn = Pti ^- (19) 

and Fx, F2 and F3 are "structure functions" describing, in the most general way 
compatible with relativistic invariance, the unknown structure of the proton. 
The antisymmetric e tensor shows that the F3 structure function is parity vi- 
olating; we will ignore it for the moment, restricting the discussion to low-Q 2 
neutral current events where the effects of W and Z exchange can be neglected. 
Contracting W^ v with the leptonic tensor used to obtain Eq. ^ leads to 

= 2+trl \ [(s + u) 2 xF 1 (x,Q 2 )-suF 2 (x,Q 2 )] , (20) 
dtdu s z t^(s + u) 1 

which can now be compared with Eq. [l7]. Substituting Eq. |l6| and evaluating the 
integral using the S function leads to 

(**\ -Vf(x) W ^ (s2 + ^ (21) 

\ d tdu) e -2?U X) sH 2 s + u ■ (21) 
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Comparison with the right-hand side of Eq. |2(j shows that, since s and u are 
continuous variables, the two equations can only be consistent if the coefficients 
of the s 2 + u 2 and su terms are equal, i.e. 

2xF 1 (x,Q 2 ) =Y^f i {x)x4 i = F 2 (x,Q 2 ). (22) 

i 

This relation between F\ and F2 is known as the Callan-Gross M relation. Equa- 
tion]^ also implies that F 2 is a function of x only, a phenomenon known as "scal- 
ing" . This was clearly observed in the original SLAC experiments, as shown in 
Fig. H and is also clearly visible in the HERA data shown in Fig. || at similar 
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Fig. 4. The vW 2 {= F2) structure function at u> = 1/x = 4 as a function of Q 2 
as measured by the SLAC-MIT group . Data taken at four different scattering 
angles are shown. All data is consistent with being independent of Q 2 . 



values of x. However, when one looks at other values of x, it is clear that scaling 
becomes progressively more and more violated. 

The phenomenon of scaling violation is one of the clearest manifestations of 
Quantum Chromodynamics and is caused by gluon radiation from the struck 
quark. This radiation is accompanied by a transfer of energy to the emitted 
gluon, which leads to a shift of the average quark x to lower values. The emitted 
gluon can also split into further quark-antiquark pairs which are also at low x. 
Thus, since the gluon bremsstrahlung depends on Q 2 , the cross section develops 
a strong Q 2 and x dependence. This can conveniently be taken into account by 
re-writing Eq. ^o| as a function of x and Q 2 using the Jacobian, (s +u)/x, as 



dxdQ 2 



2ira 2 



2na 2 



(s + u) '■ 



xF 1 (x,Q 2 )-2-F 2 (x,Q 2 



[2xy 2 F 1 + 2(1 - y)F 2 ] 



(23) 
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Fig. 5. The F-i structure function as measured by the HI and ZEUS experiments 
for bins at high function of Q 2 . The bins centred around x = 0.25 are 

where scaling was originally observed in the SLAC experiments. Clear scaling 
violation is observed in the HERA data outside this region, particularly at lower 
x. 



where we have used the identity 

p-q s + u 

V= — r = ■ (24) 

p ■ k s 

Rearranging and introducing the longitudinal structure function Fl = F<i — 2xF± 
gives 

= — w [-y 2 F L +{l + (l-y) 2 }F 2 ]. (25) 



dxdQ 2 3 

The longitudinal structure function is zero in the quark-parton model since the 
quarks have zero transverse momentum. Gluon bremsstrahlung however develops 
a non-zero pt , leading to non-zero values of Fl ■ 
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Consolidating, we can re-introduce the parity-violating xF 3 term and write 
down the most general spin-averaged form for the cross-section as 

- y 2 F L (x,Q 2 )±Y_-xF 3 (x,Q 2 )], (26) 

where the ± before xF 3 is taken as positive for electron scattering and negative 
for positron scattering, Y± are kinematic factors given by 

Y± = 1 ± (1 - y) 2 , (27) 

and 8 is the QED radiative correction. 

The F 2 structure function can be expressed, in the "DIS scheme" of renor- 
malization Q in a particularly simple way as 

F 2 {x,Q 2 ) = Y, MQ 2 )[^{x,Q 2 )+xUx,Q 2 )] (28) 

i— u,d,s : c,b 

The parton distributions qi(x, Q 2 ) and q^x, Q 2 ) refer to quarks and antiquarks 
of type i. For Q 2 -C M§, where Mz is the mass of the Z° boson, the quantities 
Ai(Q 2 ) are given by the square of the electric charge of quark or antiquark i. 
Similarly, 

xF 3 (x,Q 2 ) = Y, B i(Q 2 ) [xqi(x,Q 2 ) - xqi{x,Q 2 )] (29) 

%— u,d,s,c,b 

The full forms for the A and B terms are: 

MQ 2 ) = e 2 -2e l c e v c l v P z + (cl, 2 + c e A 2 )(c\, 2 + c\ 2 )P z 2 , (30) 
Bi(Q 2 ) = -2e lC \c\P z + Ac^J v c e A J A Pz 2 , (31) 



where 



Pz " WTmi (32) 

4 = Ji-2eiW?e w (33) 

c\ = Tl (34) 

T| = +- for i = v,u,c,t (35) 

for i — e, d, s,b (36) 



2 



4 The HERA DIS data 



4.1 The 2*2 structure function at medium and high Q 2 

The determination of the structure functions of the proton is a delicate and 
painstaking process requiring an excellent understanding of the HI and ZEUS 
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detector response. This understanding has progressed to the extent that the 
accuracy of the HERA data || [l0[ is equal to that of the fixed-target ex- 
periments in the kinematic range explored by them. Of course, the HERA data 
also extends into a much larger region of x and Q 2 and matches well onto the 
fixed-target data in the region of overlap. This has already been exhibited in 
Fig. H for the higher x bins, and is further illustrated in Fig. [| for the lower x 
bins, where the large scaling violations are even more obvious. The rapid rise of 



ZEUS+H1 




Fig. 6. The F2 structure function as measured by the HI and ZEUS experiments 
for bins at low a; as a function of Q 2 . 

the structure function at low x can be clearly seen in Fig. ^, which shows F2 in 
three Q 2 bins as a function of x. 

4.2 Next-to-leading-order QCD fits 

The precise measurement of the proton structure at low x at HERA is very 
sensitive both to the details of the evolution in QCD of the density of gluons and 
to the value of the strong coupling constant, a s , which determines the probability 
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Fig. 7. HI and ZEUS data on the structure function shown in three bins of 
Q 2 as a function of x. The steep rise of the structure function at low x is clearly 
apparent. 



of gluon emission. This sensitivity has been exploited by both ZEUS and HI. 
Each experiment has made a global QCD fit to its own data plus some or all of 
the fixed-target DIS data. There is reasonably good agreement in general terms 
between the experiments, although each experiment has a rather different fit 
procedure as well as a different choice of the fixed-target data. The quality of 
the fits in both experiments is excellent, as demonstrated by the curves shown 
in Figs. H - 0- The results for the density of the gluon are shown in Fig. ||. 

I now use the ZEUS NLOQCD fit to illustrate some points of interest. Fig- 
ure H illustrates the evolution of the gluon density and that of the sea as a 
function of Q 2 . While at medium Q 2 the sea density lies below that of the gluon 
and follows its shape, at low Q 2 it is higher than the gluon. This makes the 
normal interpretation, that the sea is driven by gluon splitting, rather difficult 
to maintain. 

Figure [lO] shows that Fl also begins to behave strangely at low Q 2 , becoming 
very flat and at the lowest values of Q 2 becoming negative, although the size of 
the uncertainties still allow it in principle to remain positive. The gluon density, 
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Fig. 8. The gluon density in the proton as measured by ZEUS (red shaded band) 
and HI jlfl] (yellow and blue shaded bands) as a function of x in three bins of Q 2 . 
The functional form used by the two collaborations in the gluon fit is somewhat 
different and is shown in the legend. 



which is directly related to Fl in QCD, certainly becomes negative, which is 
somewhat difficult to interpret in QCD. However, it is not an observable, whereas 
Fl is, so that the tendency for Fl to become negative at low Q 2 implies a break- 
down in the QCD paradigm. We will return to this discussion in Section 5.3. The 
quality of the HERA data is now so high that it alone can give constraints on 
the parton densities scarcely less good than the fits that include also fixed-target 
and other data. However, it is necessary to make 

some simplifying assumptions, particularly for the high- a; valence behaviour, 
in such fits. The results are illustrated in Fig. |HJ which shows the gluon density 
arising from such a fit, with its associated uncertainty band, in comparison with 
the result from the standard fit; the general behaviour is similar, with a somewhat 
larger uncertainty, particularly at low Q 2 . 



4.3 The determination of a s at HERA 

Another output of the NLOQCD fit is a value of a s ; the results from the two ex- 
periments are shown in Fig. |l^, labelled as "NLO-QCD fit" . The value obtained 
by ZEUS is 

a s (M 2 z ) = 0.1166 ± 0.008 ± 0.0032 ± 0.0018, 
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X 

Fig. 9. The gluon density in the proton compared to that of the quark-antiquark 
sea in bins of Q 2 as a function of x. 

where the uncertainties derive from statistical and other uncorrelated exper- 
imental uncertainties, correlated experimental uncertainties, normalisation un- 
certainties and the error related to omissions and simplifications in the NLOQCD 
model. The value HI obtain from their fit to their and the BCDMS data is 

a s {M 2 z ) = 0.1150 ± 0.0017i°:o°o5 ± 0.005, 

where the first error source takes account of all experimental uncertainties, the 
second takes account of the construction of the NLOQCD model and the final 
uncertainty results from the variation in the factorisation and renormalisation 
scale. The values of a s obtained by the two experiments are in good agreement. 

Also shown in Fig. ^ are a variety of other high-precision measurements 
of a s that can be made at HERA using a variety of techniques. These include 
classic methods such as the rate of dijet + proton-remnant production compared 
to that of single jet plus remnant, the sub jet-multiplicity evolution inside jets 
and the shape of jets. Many of these give excellent precision, comparable to the 
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Fig. 10. The F L structure function as predicted by the ZEUS NL0QCD fit as 
a function of x in bins of Q 2 . 

world average [l3| . The dominant uncertainty is usually theoretical and arises 
from the lack of predictions at next-to-next-to-leading order. 

4.4 The charm-quark structure function, F£ 

In addition to the fully inclusive structure functions discussed above, both ZEUS 
and HI can identify that fraction of F2 that arises from charm production, Fg- 
This is achieved by looking for the decay mode D* — > Dir, in which, because the 
mass difference between the D* and D is only just larger than the pion mass, the 
daughter pion has a very small momentum. The mass difference between the D* 
candidate and the D candidate can therefore be measured very accurately, allow- 
ing sufficient suppression of the combinatorial background that the charm signal 
can be cleanly identified. The structure function can then be unfolded from the 
measured differential cross section using models to correct for the unmeasured 
parts of the phase space. The measurements of the semi-inclusive charm struc- 
ture function, F£, made by both experiments [ |l4| , fl5| using this technique are 
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Fig. 11. The gluon density as a function of x in Q 2 bins from an NLOQCD fit 
using only ZEUS data. 

shown in Fig. [k| 

The data is still of rather limited statistical precision. Since the charm quark 
is produced predominantly via boson-gluon fusion, F^ is driven by the gluon 
density and thus rises steeply as x falls. 

4.5 The Fl structure function 

Since in principle both F2 and Fl are unknown functions that depend on x and 
Q 2 , the only way in which they can be separately determined is to measure the 
differential cross section at fixed x, Q 2 and at different values of y, since as shown 
in Eq. ^6], the effect of Fl is weighted by y 2 whereas F2 is weighted by 1 + (1 — y) 2 . 
However, since Q 2 — sxy, fixed x and Q 2 implies taking measurements at differ- 
ent values of s. This can certainly in principle be accomplished by reducing the 
beam energies in HERA. However, the practical difficulties for the experiments 
and the accelerator inherent in reducing either the proton or electron beam en- 
ergy, or both, by a factor sufficient to permit an accurate measurement of Fl 
mean that it has not to date been attempted. An alternative way to achieve 
the same end is to isolate those events in which the incoming lepton radiates 
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Fig. 12. Values of the strong coupling constant as determined at HERA. Each 
different measurement is displaced vertically for ease of visibility; each value 
arises from a different method as briefly indicated in the legend. The reference 
for published results is shown below the method label. The world average as 
calculated by the Particle Data Group [|l2| and by Bethke |E| are shown at the 
bottom of the figure. 

a hard photon in advance of the deep inelastic scattering, thereby reducing the 
effective collision energy. Unfortunately, the acceptance of the luminosity taggers 
typically used to detect such photons is sufficiently small and understanding the 
acceptance sufficiently difficult that no result has been obtained as yet. 

In the absence of any direct determination, the HI collaboration has used 
events at very large values of y to make an indirect measurement of Fl- The 
determinations of F2 rely on the fact that most of the measurements are made 
at values of y sufficiently small that the effects of Fl are negligible; at higher y, 
a QCD estimate of Fl, which is normally a small fraction of F2, is subtracted. 
The HI collaboration inverts this procedure by isolating kinematic regions in 
which the contribution of Fl is maximised and then subtracts off an estimate of 
F<z extrapolated from lower y. 

The method used by HI employs the derivatives of the reduced cross section 
with respect to lny. The reduced cross section can be expressed as 

xQ A d 2 a 

2ira 2 Y + ■ {I + 6) dxdQ 2 ~ ^ 
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Fig. 13. Values of the charm structure function, F 2 from the HI and ZEUS 
experiments in bins of Q 2 as a function of lnx. The blue squares show the ZEUS 
using the D* —> Kitit b decay mode, whereas the purple triangles show a ZEUS 
determination using the scmileptonic decay of the D. The shaded curves show 
the predictions from the NLO QCD fit to the inclusive F 2 data by HI. 



= F 2 {x,Q 2 )-^F L {x,Q 2 ), (37) 

which, when differentiated leads to 

do r dF 2 2y\2-y) _ f_ 3Fl_ 
dlny dlny Yf L Y+ dlny 1 ' 

which gives improved sensitivity to Fl via the stronger y dependence at the 
cost of involving derivatives of ay, F% and Fl, the quantity to be measured. It is 
instructive to consider various restrictions: 

• Small y - here da r /dhiy ~ dF^/dhiy. For low x, F 2 can be well approxi- 
mated by: 

F 2 oc x- x cx y x (39) 
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so that: 



which can be expanded as: 
dF 2 



dF 2 
dlny 



V 



d In y 



x 



Xe xlny ~ A(l + Alny...) 



(40) 



(41) 



provided Xhiy is small. From this it is clear that da r /d\ny is linear in 
lny; 

• Fl = - for all y, da r /d\ny is linear in lny for the same reason as above; 

• Fl and large y - da r /dh\y is non-linear in \ny and the deviations are 
proportional to Fl and its logarithmic derivative; 
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Fig. 14. The logarithmic derivative as a function of y in Q 2 bins. The 

curves represent the results of the HI NLO QCD fit with differing assumptions 
about Fl as shown in the legend. 

These features can be seen in the preliminary HI data of Fig. |TJ. At the 
largest values of y, the deviation from linearity implies that Fl is non-zero. 
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Although it is in principle possible to solve the differential equation for Fl 
implied by Eq. [38[ in practice the data are insufficiently precise and the QCD 
expectation is that the derivative of Fl is negligible. The uncertainty in this 
assumption is included in the systematic error. 




Fig. 15. Preliminary HI estimate of Fl. The Fl values obtained are plotted in 
Q 2 bins as a function of x. Also shown are earlier bins at higher x from the 
SLAC and NMC experiments. The curves with error bands are the predictions 
of the HI NLO QCD fit. 

The results are shown in Fig. [l^, together with earlier determinations from 
SLAC JuJ, NMC and BCDMS || ||. The curve is the result of an NLO 
QCD fit to the HI data deriving from the F% determination, i.e. by deriving the 
gluon and quark distributions from scaling violations and then calculating Fl 
using QCD. The QCD prediction is in good agreement with the HI estimate. 

5 Deep inelastic scattering at low x 

Until now we have concentrated on hard processes, in which Q 2 has been large 
and where QCD has shown itself to be applicable. However, the data, in partic- 
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ular the ZEUS "BPT" data p0[ , give access to very low Q 2 and x regimes, in 
which the strong interaction becomes very strong and perturbative QCD would 
be expected to break down. This kinematic region has traditionally been un- 
derstood in terms of Regge theory. The study of how and where the transition 
between these two regimes occurs is very interesting, not only intrinsically but 
also because of the insight it gives us into links between the apparently rather 
different processes of diffraction and deep inelastic scattering. In addition, the 
access given to very low values of x at these small Q 2 in principle gives sensitivity 
to the mechanism of QCD evolution. Given the steep rise in the parton densities 
as x falls, the data at the lowest x values may also be sensitive to high-density 
effects, such as parton recombination, sometimes known as saturation. 



5.1 QCD evolution 



We assume that the parton distribution functions, /, satisfy the schematic equa 
tion: 

df a s {n 2 ) 



(42) 



d In fi 2 2tt 

where [i represents the renormalisation scale and V is a 'splitting function' that 
describes the probability of a given parton splitting into two others. This equa- 
tion is known as the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) Equa- 
tion pi], E2, E3L B4j. There are four distinct Altarelli-Parisi (AP) splitting func- 



tions representing the 4 possible 1 



2 splittings and referred to as P qq , P gq , Pg, 



and P gg . The calculation of the splitting functions in perturbative QCD in Eq. |42| 
requires approximations, both in order of terms which can be taken into account 
as well as the most important kinematic variables. The generic form for the 
splitting functions can be shown to be [B5|: 



xV(x, a s 



oo 

£ 

n=0 



771=0 



h) I - 

X 



x¥ n \x) 



(43) 



where V (x) are the x- finite parts of the AP splitting functions and 4m' are 
numerical coefficients which can in principle be calculated for each splitting 
function. In the axial gauge, leading InQ 2 terms arise from evolution along the 
parton chain that is strongly ordered in transverse momentum, i.e. 

Q 2 > *?,„ > fc t 2 „_! > . . . 

Leading-order DGLAP evolution sums up (a s lnQ 2 )™ terms, while NLO sums 
up a s (a s lnQ ) n terms, which arise when two adjacent transverse momenta 
become comparable, losing a factor of lnQ 2 . 

In some kinematic regions, and in particular at low x, it must become essen- 
tial to sum leading terms in In 1 /x independent of the value of In Q 2 . This is done 
by the Balitsky-Fadin-Kuraev-Lipatov ||| §7|, |j| |§ (BFKL) equation, which 
governs the evolution in x at fixed Q 2 . The leading-order terms in (a s lnl/a;) n 
arise from strong ordering in x, i.e. 



x <C x n <C x n -i <C 
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One of the most important goals of HERA physics is the search for experimental 
effects that can be unambiguously attributed to BFKL evolution. 

Generally, however, QCD coherence implies angular ordering. To see the im- 
plications of this it is more convenient to work with unintegrated parton density 
functions, f(x, fc 2 , /z 2 ), where jj, is the scale of the probe. There are now two 
hard scales, kt and complicated QCD evolution, known as the Ciafaloni, Catani, 
Fiorani and Marchesini |jj ^ §| (CCFM) evolution equation. The DGLAP 
and BFKL equations can then be seen to be two limits of angular ordering. In 
the DGLAP collinear approximation, the branching angle, 9, where 9 ~ kt/ki, 
grows since k t grows; while for BFKL evolution, 6 grows because ki oc x falls. 

Figure |l6| shows the lnl/x - InQ 2 plane at HERA, together with schematic 
indications of the directions in which GLAP, BFKL and CCFM evolution is 
expected to be most applicable. Also indicated on the figure are schematic in- 




Fig. 16. Schematic diagram showing different regions of the lnl/x and InQ 2 
plane and the evolution equations expected to hold therein. The line marked 
'saturation' represents the boundary between GLAP evolution and evolution 
governed by the GLR equation. The 'size' of partons is also indicated in differing 
kinematic regions. 

dications of both the 'size' and density of partons in the proton in different 
kinematic regions. 

The transverse size of the partons that can be resolved by a probe with 
virtuality Q 2 is proportional to 1/Q, so that the area of the partonic 'dots' in 
Fig. H falls as Q 2 rises. For particular combinations of parton size and density, 
the proton will eventually become 'black' to probes, or, equivalently, the compo- 
nent gluons will become so dense that they will begin to recombine. The dotted 
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line labelled 'Critical line - GLR' refers to the boundary beyond which it is ex- 
pected that such parton saturation effects will become important, i.e. the region 
in which partons become so densely crowded that interactions between them 
reduce the growth in parton density predicted by the linear GLAP and BFKL 
evolution equations. The parton evolution in this region can be described by the 
Gribov-Levin-Riskin p3| , p4[ equation, which explicitly takes into account an 
absorptive term in the gluon evolution equation. Naively, it can be assumed |35| 
that the gluons inside the proton each occupies on average a transverse area of 
irQ~ 2 , so that the total transverse area occupied by gluons is proportional to 
the number density multiplied by this area, i.e. TrQ~ 2 xg(x, Q 2 ). Since the gluon 
density increases quickly as x falls, and the gluon 'size' increases as in the 
region in which both x and Q 2 are small, saturation effects ought to become 
important. This should occur when the size occupied by the partons becomes 
similar to the size of the proton: 

xg(x,Q 2 )^ = nR 2 (44) 

where R is the radius of the proton, (~ 1 fm ~ 5 GeV -1 ). The measured values 
of xg(x, Q 2 ) imply that saturation ought to be observable at HERA [[36| at low 
x and Q 2 , although the values of Q 2 which satisfy Eq. |44| are sufficiently small 
that possible non-perturbative and higher-twist effects certainly complicate the 
situation. Of course, 

it is also possible that the assumption of homogenous gluon density is in- 
correct; for example, the gluon density may be larger in the close vicinity of 
the valence quarks, giving rise to so-called 'hot spots' j37|, which could lead to 
saturation being observable at smaller distances and thereby larger Q 2 . 

5.2 Interpretation and Models 

The region of low Q 2 and low x is one in which perturbative QCD meets and 
competes with a large variety 

of other approaches, some based on QCD, others either on older paradigms 
such as Regge theory or essentially ad-hoc phenomenological models. In a pre- 
vious article Q I gave a quite detailed review of these models, and given that 
there have not been major developments here, I refer the reader to that for 
further information. 

5.3 QCD fits 

The extension of the kinematic range and the high-precision data on F% from 
HERA provided a substantial impetus to the determination of parton distri- 
bution functions via global fits to a wide variety of data. The major current 
approaches are due to the CTEQ group §f] and Martin et al. (MRST) Q. 
In general both groups fit to data from fixed target muon and neutrino deep 
inelastic scattering data, the HERA DIS data from HERMES, HI and ZEUS, 
the W-asymmetry data from the Tevatron as well as to selected process vary- 
ing from group to group such as prompt photon data from Fermilab as well as 
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high- Et jet production at the Tevatron. The different data sets give different 
sensitivity to the proton distributions, depending on the kinematic range, but 
together constrain them across almost the whole kinematic plane, with the pos- 
sible exception of the very largest values of x, where significant uncertainties still 
remain pj[ . 

The approaches of CTEQ and MRST are basically similar, although they dif- 
fer both in the data sets used as well as in the fitting procedure and the technical 
details of the theoretical tools used, e.g. the treatment of heavy quarks in DIS. 
In their latest fits, CTEQ prefer to omit the prompt photon data because of the 
uncertainties in scale dependence and the appropriate value for the intrinsic kx 
required to fit the data. Instead they use single-jet inclusive Ej~ distributions to 
constrain the gluon distribution at large x. In contrast, until their most recent 
publications, MRST retained the prompt photon data, giving alternative PDFs 
depending on the value for the prompt-photon intrinsic kr used. Both groups 
parameterise the parton distributions in terms of powers of x and (1 — x) lead- 
ing to fits with many free parameters. The MRST NLO parameterisation of the 
gluon is shown below as an example: 



xg = A g x x ° (1 - x) ri! > (1 + e g y/x + j g x) (45) 

where A g ,X g ,r]g, e g and j g are free parameters in the fit. 

The seminal work of Botje in producing PDFs with associated error matrices 
for the first time p2] has led to similar fits being produced by other groups. 



CTEQ 1 43 have produced a global fit with associated errors, while both HI 



and ZEUS |45[ have produced their own fits using DIS data only, as discussed 



in Section 4.2 



These fits allow one to see very graphically the salient features of the QCD 
evolution that have been discussed above. For example, the ZEUS NLOQCD fit 
shown in Fig. |9|has already been discussed in Section 4.2. Here we saw that the 



strange behaviour both of the gluon vs. sea densities, as well as Fl, as a function 
of Q 2 showed that the normal QCD interpretation may well be breaking down 
at Q 2 around 1 GeV 2 . 

One possible reason for the problems with the interpretation of the NLO QCD 
fits could be the need for higher-order QCD fits i.e. next-to-next-to-leading-order 
(NNLO) fits. The first steps in implementing such fits have already begun; some 
moments of the NNLO splitting functions have been calculated fiJJ . Using this 
with other available information, van Neerven and Vogt [^7], Q have produced 
analytical expressions for the splitting functions which represent the slowest and 
fastest evolution consistent with the currently available information. The MRST 
group has recently used this information to investigate NNLO fits to the available 
data |49f . Such an analysis requires some changes to the parameterisations used, 
so that for example the NLO parameterisation of the gluon of Eq. |4^ becomes: 

xg{x, Ql) = A g x~ x ° (1 - xf' (1 + e g Vx~ + l 9 x) - A'g x~ w » (1 - x)< , (46) 

primarily in order to facilitate a negative gluon density at low x and low Q 2 , 
which, as we have seen, is preferred by the fits, even at NLO. The results of the 
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Fig. 17. The MRST 'central' NNLO fit to DIS data. The solid line shows the 
NNLO fit, while the NLO fit is shown by the dashed line and the LO fit by the 
dotted line. The data are from HI, ZEUS and the fixed target experiments and 
are plotted in x bins as a function of Q 2 with an additive constant added to the 
data of each x bin to improve visibility. 



'central' fit, between the extremes of the van Neerven-Vogt parameterisation, is 
shown in Fig. [H]. 

There are also changes of the LO and NLO fits with respect to earlier publi- 
cations, in as much as MRST now follow CTEQ in using the Tevatron high-£V 
data rather than the prompt-photon data, and HERA F-i data has been included 
in the fit. There is a marked improvement in the quality of the fit in the pro- 
gression LO — > NLO — > NNLO, in particular in terms of the NMC data. The 
size of higher-twist contributions at low x also decreases, so that at NNLO is it 
essentially negligible. The effect of going to NNLO on the PDFs themselves is 
highly non-trivial. This is illustrated in Fig. [l^, where the quite major changes 
in Fl, particularly at low x, are evident. There is also a large variation depend- 
ing on the choices made in the parameter space allowed by the partial NNLO 
ansatz. Indeed, the GLAP approach is not convergent for Q 2 < 5 GeV 2 , which 
may well be due to the neglect of important lnl/x contributions. However, the 
instability seen at low Q 2 soon vanishes at higher Q 2 . 
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Fig. 18. The Fl structure function from the MRST fits, taking into account 
part of the NNLO corrections in four bins of Q 2 as a function of x. The solid 
line shows the 'average' of the parameter space available to choose the NNLO 
parameters, while the dashed-dotted lines show the two extreme possibilities. 
The NLO fit is indicated by the dashed line while the LO fit is indicated by the 
dotted line. 



Thorne has investigated the question of incorporating In 1/x terms in the 
splitting functions by incorporating the solution of the NLO BFKL kernel using 
a running coupling constant |^0|, The inclusion of the BFKL terms does 
indeed give an improved fit compared to the 'central' NNLO fit, particularly at 
the lowest Q 2 and x. This may be one of the first unambiguous indications of 
the importance of BFKL evolution. 

5.4 Data in the transition region 

The approximate position of the transition between data that can be described 
by pertrubative QCD evolution and that which require the Regge approach can 
be seen in Fig. |l9|. 

For Q 2 > 1 GeV 2 , the data are roughly independent of Q 2 , whereas at lower 
Q 2 they fall rapidly, approaching the Q~ 2 dependence that would be expected 
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Fig. 19. ZEUS BPT data on F2 in bins of y as a function of Q 2 . Also shown are 
earlier ZEUS data as well as data from HI. The solid line shows the results of 
the ZEUS Regge fit to the form of Eq. [It], while the dotted line shows the result 
of the ZEUS NLO QCD fit. 



in the limit Q 2 — > from conservation of the electromagnetic current. Although 
QCD gives a good fit to the data down to Q 2 ~ 1 GeV 2 , below that it is necessary 
to use a Regge-based fit of the form 

^-(&)(^)(-(?r^(?r). 

(47) 

where Ar,, Ajp and Mq are constants and «r and ajp are the Reggeon and 
Pomeron intercepts, respectively. Regge theory is expected to apply at asymp- 
totic energies. The appropriate energy here is W, the centre-of-mass energy of 
the virtual photon-proton system, given by Eq. || Since, at low x, W 2 ~ 1/x, it 
would be expected that Regge fits would be applicable at very low x and Q 2 . 

The complete ZEUS data over six orders of magnitude in x and Q 2 are shown 
in x bins as a function of InQ 2 in Fig. together with fixed target data from 
NMC and E665, which extends the range in the direction of medium x and Q 2 . 

The availability of this very wide range of precise data makes possible qual- 
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Fig. 20. Compilation of ZEUS Fi data in x bins as a function of Q 2 . Each x 
bin is shifted by an additive constant for ease of visibility. Data from NMC and 
E665 are also shown. The dotted lines show lines of constant W, while the solid 
lines are fits to the form of Eq. 

itatively new investigations of models that describe Fi- Since the logarithmic 
derivative of F 2 is directly proportional to the gluon density in leading-order 
QCD, which in turn is the dominant parton density at small x, its behaviour 
as a function of both x and Q 2 is important. The solid curves on the figure 
correspond to fits to a polynomial in In Q 2 of the form 

F 2 = A(x) + B(x) (log 10 Q 2 ) + C(x) (log 10 Q 2 f , (48) 

which gives a good fit to the data through the entire kinematic range. The dotted 
lines on Fig. ^ are lines of constant W. The curious 'bulging' shape of these 
contours in the small- a; region immediately implies that something interesting is 
going on there. Indeed, simple inspection of Fig. ^ shows that the slope of F 2 
at constant W begins flat in the scaling region, increases markedly as the gluon 
grows and drives the evolution of F2 and then flattens off again at the lowest x. 
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Figure ^T] shows the logarithmic derivative evaluated at (x, Q 2 ) points along 
the contours of fixed W shown on Fig. according to the derivative of Eq. fl8| , 
viz.: 

2 — = B(x)+2C(x)\og w Q 2 , (49) 

° iO 610 V 

where the data are plotted separately as functions of InQ 2 and In a;. The turn- 
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Fig. 21. The logarithmic derivative of the ZEUS Fi data in six bins of W, plotted 
as a function of Q 2 and x. 

over in the derivatives in all W bins is marked. Within the framework of pQCD, 
the interpretation of such an effect is that the growth of the gluon density at 
low x is tamed as Q 2 and x fall. This behaviour can also be seen more clearly 
in Fig. which shows a three-dimensional plot of the derivatives as a function 
of both In a: and InQ 2 , obtained from a parameterisation of the DIS data. 

Such an effect is by no means necessarily an indication of deviations from 
the standard DGLAP evolution. It can be seen from Fig. ^2] that, for example, 
no turnover effect occurs in bins of constant Q 2 . Nevertheless, the features of 
Figs. |2l| and ^2] can be explained as a natural consequence of parton saturation 
or shadowing. These effects can be naturally discussed in "dipole models" Q, 
which often explicitly take into account parton-saturation effects. In such models, 
the "standard" picture of deep inelastic scattering in the infinite-momentum 
frame of the proton is replaced by an equivalent picture produced by a Lorentz 
boost into the proton rest frame. In this frame, the virtual photon undergoes 
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Fig. 22. The logarithmic derivative of the ZEUS F% data plotted as a surface in 
three dimensions versus both Q 2 and x. The curves show lines of constant W . 

time dilation and develops structure far upstream of the interaction with 
the proton. The dominant configurations of this structure are qq and qqg Fock 
states, which interact with the proton as a colour dipole. The higher the Q 2 of 
the interaction, the smaller the transverse size of the dipole. For small x, the deep 
inelastic process can be considered semi-classically as the coherent interaction 
of the dipole with the stationary colour field of the proton a long time after 
the formation of the dipole. As an example, the model of Golec-Biernat and 
Wusthoff (GBW) |52|, [33] is shown in Fig. |2^, together with the results of the 
ZEUS NLOQCD fit. 

It can be seen that the GBW model reproduces the basic features of the log- 
arithmic derivative plot reasonably well. However, so does the ZEUS NLOQCD 
fit, so that no firm conclusion can be drawn on the existence or otherwise of 
saturation effects. 

+ + 
6 Diffraction 

Diffractive DIS is that subset characterised by a hard interaction between the 
proton and the exchanged virtual photon that nevertheless leaves the proton 
intact. Such interactions are normally thought of in Regge theory as proceeding 
via the exchange of a colourless particle with the quantum numbers of the vac- 
uum, known as the Pomeron. The Pomeron under some circumstances can be 
considered to develop its own partonic structure, analogous to the proton, which 
can be parameterised using diffractive DIS data. 

One of the most attractive features of dipole models such as that of GBW 
discussed in the previous section is the rather natural way in which they can lead 
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Fig. 23. The logarithmic derivatives of the ZEUS F% data as shown in Fig. 
The dotted curves show the predictions of the BGW model for the five highest W 
bins; since these predictions do not include QCD evolution effects, they should 
become progressively less accurate as Q 2 increases. The solid curves show the 
predictions of the ZEUS NLO QCD fit for all W bins; they are only shown down 
to Q 2 ~ 3 GeV 2 . 

to a unified description of diffraction and deep inelastic scattering. The QCD in- 
terpretation of the Pomeron is that it is equivalent to the exchange of two gluons 
in a colour-singlet state. In this picture, therefore, diffraction can be considered 
as a subset of fully inclusive DIS, which sums over all possible exchanges between 
the dipole and the proton, dominantly one- and two-gluon exchange in a colour 
octet, in contrast to the colour-singlet exchange that dominates inclusive DIS. 
This deep connection between these two processes leads to non-trivial predic- 
tions which do indeed seem to be at least qualitatively in agreement with the 
data. This is illustrated in Fig. ^J. This figure is surprising for several reasons. 
It demonstrates that the diffractive cross section has the same W dependence as 
the total cross section. To the extent to which the diffractive cross section can 
be related to the elastic cross section, one would have expected from the Optical 
Theorem that the ratio would have a power-law dependence on W, as indeed 
would also be expected from Regge theory via the exchange of a Pomeron. A 
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Fig. 24. The ratio of the diffractive to total cross section in four Q 2 bins as 
a function of W. The curves show the predictions from the Golec-Biernat & 
Wusthoff model. 

strong W (~ 1/x) dependence is also expected in QCD models, since the to- 
tal cross section is dominated by single-gluon exchange, whereas diffraction is 
dominated by two-gluon exchange. The other surprise is the fact that the GBW 
model gives a rather good qualitative representation of the data. 

ZEUS has also investigated the behaviour of this ratio as Q 2 — > |54| . Fig- 
ure ^5] shows the diffractive structure function, x-pF^^ (the analogue to F%, 
integrated over t) multiplied by , the fraction of the protons momentum car- 
ried by the Pomeron, as a function of Q 2 in various bins of W and in two bins 

D(3 

of Mx, the mass of the hadronic system other than the proton. The F 2 y data 
points are determined using two methods, one of which requires the observation 
of a large rapidity gap in the proton-beam direction, while the other uses the 
ZEUS Leading Proton Spectrometer (LPS) [p5| . This device is an array of six 
stations of silicon-strip detectors placed in Roman pots downstream of the in- 
teraction point in the proton beam direction. It uses the HERA beam elements 
to form a magnetic spectrometer to analyse the leading proton from diffractive 
interactions. Although its acceptance is of necessity small, it avoids the low-mass 
proton dissociative background endemic with other forms of identifying diffrac- 
tive interactions. The similarity of Fig. |25| to Fig. |l9| is striking, again illustrating 
the Q~ 2 falloff enforced by electromagnetic current conservation. 
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Fig. 25. xjpF 2 in bins of W and Mx as a function of Q 2 . The curves refer 
to the model of Bartels et al. |36) 

6.1 Vector meson production 

The exclusive production of vector mesons is both a very simple laboratory to 
study many aspects of diffraction as well as a process in which dipole models 
are likely to be particularly appropriate. Figure shows the photoproduction 
cross sections for a variety of vector mesons as a function of W, as well as the 
total cross section. The relatively slow rise of the total cross section with Q 2 is 
indicative of the dominance of soft processes. The J/ip cross section clearly has 
a much steeper rise with Q 2 . This is more clearly seen in Fig. |2?], where the cross 
section for several different Q 2 is shown as a function of W. It would appear that 
the mass of the J/ ip is sufficiently large that it gives rise to a hard scale even at 
Q 2 ~ 0. 

In contrast, the mass of the p is small and quite large values of Q 2 need 
to be reached before the W dependence rises to the values associated with hard 
processes. This can be seen in Fig. |2^, which shows the values of the fit to a power 
law in W as a function of Q 2 . An appropriately hard scale seems to pertain for 
Q 2 > 5. 

This behaviour, in which either Q 2 or the mass can act as a hard scale, leads 
to the obvious question of whether a combination of these two quantities can 
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Fig. 26. The total photc-production cross section together with the photopro- 
duction cross sections for a variety of vector mesons as a function of W from 
HI, ZEUS and several fixed-target experiments. The dotted lines show W to the 
indicated powers. 

also give a hard scale. Figure ^8] shows the data on /?, tfi and J/ip production 
from HI and ZEUS plotted against Q 2 + M 2 . There is indeed a tendency for 
the data for all the vector mesons to lie on a universal curve. However, other 
more detailed comparisons |57j ] show that some differences do remain between 
the different species even when plotted against Q 2 + M 2 . 

To the extent that we can model diffraction by the exchange of a colourless 
two-gluon state, we would expect a difference between the W dependence of 
vector meson production and inclusive DIS. This is illustrated in Fig. ^ shows 
the values of fits to the JU-dependence of the inclusive DIS and vector-meson 
cross sections against Q 2 and Q 2 + M 2 , respectively. For Q 2 + M 2 greater than 
about 5 GeV 2 , the value of 5 is indeed about twice that at the same value of 
Q 2 in inclusive DIS, as would be expected in the simple picture of two-gluon 
exchange. 

Provided that we have a scale sufficiently hard that pQCD is applicable, 
we can also use vector meson production to probe the gluon density in the 
proton. Since the cross section is proportional to the gluon density squared, 
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Fig. 27. The ZEUS and HI J/tp cross section for several bins of Q 2 as a function 
of W. The dotted lines show the predictions of the model of Frankfurt, Koepf 
and Strikman , the dashed lines show the predictions of Martin, Ryskin and 
Teubner |55| and the solid lines show fits to the form of W s . The inset shows 
the S obtained from these fits as a function of Q 2 . 



such a determination should in principle be much more sensitive to the gluon 
density than, for example, scaling violations in deep inelastic scattering. This is 
illustrated in Fig. [30|, which shows ZEUS and HI data for J/ip photoproduction 
as a function of W. We have seen already that the J/ip mass is sufficiently large 
to guarantee that pQCD is applicable even in photoproduction. The quality of 
the data is sufficiently high that it is in principle sensitive to the gluon density. 
Unfortunately, the wave-function of the J/tp must be modelled, which leads to 
substantial uncertainty in the model predictions, so that this data has not as 
yet been used in global fits to constrain the gluon density. The good agreement 
between at least some version of the models and the data, as shown in Fig. 
does however indicate that the gluon determined in the global fits to DIS data is 
indeed also able to explain the dynamics of this completely different diffractive 
process. 

Finally, it is interesting to investigate whether t can also provide a hard scale 
for pQCD. ZEUS has precise data out to —t ~ 12 GeV 2 for p production and 
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Fig. 28. a)Results of a fit to the ZEUS p data of the form W s . The value of the 
exponent is plotted as a function of Q 2 . b) The exponent obtained from fits of 
the form W s to the vector meson data from ZEUS and HI. The value of the 
exponent S is plotted against Q 2 + AI 2 . 



also to beyond 6 GeV 2 for 4> and J/ip- The ratio of the cross sections for and 
p starts somewhat below the SU(4:) expectation but reaches it quite quickly, by 
— t ~ 3 — 4 GeV 2 . In contrast, the J/ip ratio remains below the SU(4) prediction 
for much longer, hardly reaching it even for — t ~ 6 GeV 2 . Figure |3l] shows 
the ZEUS data for all three vector mesons as a function of — t compared to 
a two-gluon exchange model and to a model in which a BFKL gluon ladder 
is exchanged It can be seen that the two-gluon model completely fails to 
reproduce the data both in magnitude and in shape, whereas the BFKL model, 
which has been fit to the ZEUS data, gives an excellent fit for all three mesons. 
Although some of the model assumptions, such as a fixed a s and a (5-function 
wave-function for the light vector mesons, are somewhat questionable, the fact 
that this BFKL model fits the data whereas conventional pQCD fails is very 
striking. It can thus be concluded that t does form an appropriate hard scale 
for pQCD calculations. However, the ratios of ZEUS vector-meson cross sections 
behave differently as a function of t compared to Q 2 , so that either these ratios 
are different at asymptotic values of the two variables, or t and Q 2 are not 
equivalent hard scales. 

6.2 Deeply virtual Compton scattering 

Deeply virtual Compton scattering (DVCS) is an interesting process in that it 
is the simplest possible non-elastic diffractive process. It consists of diffractive 
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Fig. 29. a) Values of the exponent in fits of the form W to the W dependence 
of the inclusive DIS total cross section, as a function of Q 2 . b) Values of the 
exponent in fits of the form W 5 to the W dependence of the cross section for a 
variety of vector mesons, as a function of Q 2 + M 2 . 



scattering of the virtual photon from the proton, putting the virtual photon 
onto mass shell so that the final state consists of the initial proton and positron 
plus a photon. As such, no fragmentation or hadronic wave-functions complicate 
the process. The Bethe-Heitler QED process also leads to the same final state, 
which is both a blessing and a curse. The curse is that is necessary to separate 
the two processes; the blessing is that the fact of identical final states leads to 
interference. Since the amplitude and phase of the QED process is calculable, 
in principle this opens the door to the determination of the unknown QCD 
amplitude that governs DVCS. However, we are still some considerable distance 
from this goal. 

Given that a virtual photon participates in the collision with the parton 
in the proton and that it emerges on mass shell, it is clear that the scattered 
parton that must be re-integrated into the final-state proton must undergo a 
change in its four-momentum corresponding to a change in its x value. Thus the 
DVCS process is sensitive to the so-called "skewed parton distributions inside 
the proton, which can be thought of as the cross-correlation function between 
partons of fractional momentum X\ and x%. Thus much unique information can 
be obtained by studying these processes. 

The DVCS process was first seen at HERA by ZEUS |6l). The HI collab- 
oration has published results on the observation of this process and the first 
measurement of the cross section 62 1. The background from the Bethe-Heitler 
process can be subtracted by utilising the kinematic characteristics of the two 
processes. In DVCS, the photon is normally produced at a large angle to the 
incident beam directions and the positron at a small angle, whereas, for the 
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Fig. 30. a) The ZEUS and HI data on J/tp photoproduction, together with the 
data from two fixed target experiments as a function of W. The curve shows 
a fit to the form W s . b) The ZEUS data together with curves showing the 
predictions using the stated gluon pdfs from the models of Frankfurt, McDermott 
and Strikman |6f| and Martin, Ryskin and Teubner p9[ . 



Bcthe-Heitler process, the reverse is the case. The sum of Monte Carlo gen- 
erators describing the two processes gives a good description of the kinematic 
quantities of the data, giving confidence that the DVCS cross section can be 



extracted. The HI and ZEUS cross sections are shown in Fig. 32. It can be seen 
that the cross sections are completely dominated by statistical errors and that 
both collaborations find good agreement with theoretical models, notably that 
of Frankfurt, Freund and Strikman Q and Donnachie and Dosch [Q. As more 
data is collected, the possibility of using this process to determine skewed parton 
distributions will make it a fruitful area of study at HERA II. 

There are many more areas of diffraction, in particular the diffractive struc- 
ture functions, which I have not been able to cover. The study of diffraction has 
turned out to be one of the most exciting and rich areas of HERA physics. It 
will continue at HERA II and will hopefully allow us to make progress in our 
understanding of soft interactions and the related, and fundamental, question of 
understanding confinement in QCD. 



7 Heavy quark production 

We have already touched upon the production of charm quarks at HERA in 



the discussion of the F$ in Section 4.4. In fact, copious amounts of charm are 



produced in photoproduction at HERA, so that for many channels, HERA is 
a "charm factory. The ZEUS collaboration has made several contributions to 
charm spectroscopy results: as an example, Fig. p3 shows the D*ir spectrum 
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Fig. 31. The ZEUS data on vector- meson production at high —t. The curve 
shows a fit to a BFKL-model by Forshaw and Poludniowski. 



from 110 pb" 1 of data. A rather complex set of structures can be observed 
between 2.4 and 2.5 GeV, some of which correspond to known resonances and 
some of which do not. Work continues to understand this complex area; the 
sensitivity of ZEUS in this sort of investigation is similar to that of the LEP 
experiments and CLEO. 

Both HI and ZEUS have observed beauty production in photoproduction; 
HI has also published a cross section in DIS. The identification of the b signal 
is based on the use of high-transverse-momentum leptons for ZEUS; HI also use 
tracks with large impact parameter as measured in their silicon vertex detector. 
The results as a function of Q 2 are shown in Fig. |34|. It can be seen that the QCD 
predictions are substantially below the data for all Q 2 , thus joining a pattern 
also seen in proton-antiproton and photon-photon collisions. The tendency for 
pQCD to fail to predict B cross sections, which naively would be thought to be 
an area in which it should work well, is becoming increasingly interesting. The 
advent of HERA II will make an enormous difference to the precision of this type 
of measurement and should allow a stringent test of the theory of heavy-quark 
production. 

8 High-Q 2 phenomena 

HERA provides an unique opportunity to study the electroweak interaction at 
Q 2 sufficiently high that the charged and neutral currents are of similar strength. 
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Fig. 32. a) The cross section for deeply virtual Compton scattering as a function 
of Q 2 as measured by the HI collaboration. The dotted area corresponds to 
the prediction of Frankfurt, Freund and Strikman, while the solid area is the 
prediction of Donnachie and Dosch. b) The DVCS cross section as measured 
by ZEUS. The curve shows the prediction for the elastic process by Frankfurt, 
Freund and Strikman. 



Figure 35 shows the differential cross-sections for the charged and neutral cur- 
rents as a function of Q 2 from HI and ZEUS. It can be seen that, for e~p 
interactions, these two processes become of equal strength at Q 2 ~ M§ ~ 10 4 
GeV 2 . For e + p interactions, the charged current cross-section approaches the 
neutral current cross-section, but remains below it. The features of this plot can 
be explained by inspection of Eq. together with Eqs. ^30| and [5l] below: 



d 2 c 



dxdQ 2 



cc 



G 2 F 
27 



+ Q 2 



2x{u(x) + c(x) + (1 - y) 2 (d(x) + s(x))} 
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2x{u(x) + c(x) + (1 - y) 2 (d{x) + s(x))} (51) 

For the charged current case, the smaller size of the e + p cross-section com- 
pared to e~p is related to the fact that, at high Q 2 , Eq. || implies that both 
x, y — ► 1. There are two main contributory factors to the cross-section difference 
that flow from this. First, there are twice as many u valence quarks inside the 
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Fig. 33. The D*tt spectrum from ZEUS, where the D* decays via Dw and the 
D via Kir. To improve the accuracy of the reconstruction, the effective mass of 
the Kiririr system minus the effective mass of the Kttt: system that forms the 
D* plus the D* mass is plotted. The curves show a polynomial background plus 
Gaussians for the known resonance at 2.420 GeV and the D\ at 2.460 GeV. 
Another structure at just above 2.4 GeV also seems to be visible. 



proton that can couple to W~ as d quarks that can couple to W + . Secondly, the 
(1 — y) 2 terms in Eqs. |5(] and ^l], which arise from the V — A helicity structure 
of the charged weak current, imply that the valence-quark contribution, which 
is dominant at high Q 2 , is suppressed for the positron case but not for electrons. 

The difference between the electron and positron neutral current cross sec- 
tions shown in Eq. ^6] allows the determination of the parity- violating structure 
function xF% by taking the difference of the cross sections. The results |Q are 
shown in Fig. |3^. Since its determination requires the subtraction of two quan- 
tities that are almost equal, it is dominated by statistical uncertainties, which 
are in turn dominated by the fact that the electron data sample that has so far 
been obtained at HERA is much smaller than that for positrons. 

The high-Q 2 regime is also interesting since possible new states from electron- 
quark fusion (e.g. leptoquarks) have masses given by M 2 ~ sx and since the 
sensitivity to the effects of new currents is maximised. An example of the sensi- 
tivity that can be obtained at HERA is shown in Fig. ^57|, which shows the mass 
against coupling limits for two varieties of scalar leptoquark. Both HI and ZEUS 
have comparable limits for a whole range of such states with differing quantum 
numbers. It can be seen from Fig. and it is generally the case, that for some 
states, in particular in i?-parity-violating supersymmetry models or leptoquarks, 
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Fig. 34. The ratio of the measured cross section and the theoretical prediction, 
plotted against Q 2 , from HI on beauty production in photoproduction and DIS 
and from ZEUS in photoproduction,. The band shows the uncertainty on the 
predictions of perturbative QCD. 



HERA has higher sensitivity than either LEP or the Tevatron. 

Limits on excited leptons and quarks have also been obtained by ZEUS and 
HI, which extend the limits from the LEP experiments considerably beyond the 
LEP II centre-of-mass energy. 

As well as stringent limits on new phenomena, the HERA data also show 
intriguing features which may be signatures for new physics. The HI collabo- 
ration has observed a class of events that have isolated charged leptons with 
large missing transverse momentum. Figure |3^ shows the distribution of the 
transverse momentum of the hadronic system, p?£ , against its transverse mass, 
separately for electrons (or positrons) and muons in such events. Also shown are 
the expectations from the Standard Model background, which is dominated by 
single W production. 

It can be seen that the distribution of the events is rather different to 
the Standard Model expectation. Furthermore, for the transverse mass of the 
hadronic system greater than 25 GeV, HI sees four electron and six muon events, 
compared to Standard Model expectations of 1.3 and 1.5 events, respectively. 
Unfortunately, this exciting observation is not confirmed by ZEUS, which, for 
the same cut in p?f , sees one event in each category compared to the Standard 
Model expectation of 1.1 and 1.3, respectively. Intensive discussions between the 
two experiments have not revealed any reason why HI might artificially produce 
such an excess nor why ZEUS should not observe it. It would therefore seem that 
there must be an unlikely fluctuation: cither the HI observation is an upward 
fluctuation from the Standard Model, or ZEUS has suffered a downward fluctu- 
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Fig. 35. Charged and neutral current differential cross sections for e scattering 
as a function of Q 2 from HI and ZEUS. 

ation from a signal for new physics. More data from HERA II will be required 
to resolve this puzzle. 

One possible source of an excess of events with isolated leptons with missing 
transverse momentum would be from a flavour-changing neutral current process 
producing single top quarks. Both HI and ZEUS have used the samples described 
above to put limits on the FCNC couplings of the 7 to light quark-top quark 
vertices. The results are shown in Fig. [?9|. Also shown are the limits from LEP 
and CDF, which are complementary to those from HERA, in the sense that, 
since the Z-exchange cross section at HERA is so much smaller than that for 7 
exchange, the HERA data limit only the photon coupling. 

9 HERA II physics 

Since many of the physics results discussed above, particularly those at high Q 2 , 
are statistics limited, there is a clear physics case for a significant increase in 
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Fig. 36. The xF% structure function as determined by HI and ZEUS as a function 
of x in six bins of Q 2 . 



integrated luminosity for HI and ZEUS. There are also other interesting physics 
investigations possible at HERA that have not yet been carried out. For exam- 
ple, there is a natural build-up of transverse polarisation of the lepton beam in 
HERA that occurs through the Sokholov-Ternov effect |6^] . As very successfully 
demonstrated at HERMES using gas targets, this transverse polarisation can be 
rotated into the longitudinal direction and utilised to do physics. The installation 
of spin rotators in HI and ZEUS would allow polarisation studies to be carried 
out at very much higher Q 2 . This is particularly interesting to study the chiral 
properties of the electroweak interaction. For these and several other reasons, it 
was decided to embark on a major upgrade of both the HERA accelerator and 
the HI and ZEUS detectors. The aim of the HERA II programme is to produce 
a factor of approximately five increase in luminosity and accumulate 1 fb _1 of 
data with both electron and positron collisions in both longitudinal polarisation 
states. 

The changes to the HERA accelerator include the replacement of 480 meters 
of the vacuum system and the design and installation of almost 80 magnets in 
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Fig. 37. Limits on coupling strength A versus mass Mlq for leptoquarks. The 
top plot shows limits for fermion number = leptoquarks decaying into the eq 
final state from ZEUS. The lower plot shows limits from HI for fermion number 
= 2 leptoquarks decaying into both eq and vq final states. Also shown are limits 
obtained from the Tevatron (yellow shaded area) and LEP (blue striped area). 
These leptoquark species have identical quantum numbers to squarks that violate 
R-parity. 

the region around the HI and ZEUS interaction points. In particular, supercon- 
ducting quadrupole focussing elements were inserted inside both detectors to 
reduce the beam emittance and spin rotators were installed on either side of the 
HI and ZEUS interaction regions. 

Both the ZEUS and HI detectors have undergone a massive programme 
of consolidation and repair work, as well as major detector upgrades. As an 
example, I discuss briefly the changes made to ZEUS; the general thrust of the 
upgrade is similar in the two detectors, although the details are different. 

9.1 Upgrades to ZEUS for HERA II 

The ZEUS upgrades have concentrated in three main areas: the vertex region; 
the forward direction; and the luminosity monitoring. 

9.1.1 The vertex region The tagging of the large flux of heavy quarks 
(charm and beauty) produced at HERA II can be greatly enhanced by the in- 
stallation of a high-precision charged-particle detector as close as possible to a 
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Fig. 38. Distribution of transverse mass versus the pr of the hadronic system 
for HI events containing a) isolated electrons and b) isolated muons. The dots 
show the distribution of Standard Model W Monte Carlo events corresponding 
to a luminosity 500 times that of the data. 



thin beampipe. The ZEUS MVD ]6^, |6J| consists of 20 fim pitch n-type silicon- 
strip detectors with p + -type implants. The readout pitch is 120 //m, leading to 
more than 200,000 readout channels, which are digitised by a custom-built clock, 
control and ADC system. The detectors are organised in two main groups: a "bar- 
rel", which surrounds the elliptical 2 mm-thick 1.1% of a radiation length) 
aluminium-beryllium beam-pipe; and four "wheels" , consisting of wedge-shaped 
detectors mounted perpendicular to the beam-line in the forward direction from 
the interaction point. Figure flO] shows one half of the MVD before installation 
at DESY. In the barrel region, the ladders, each of which consists of five silicon 
detectors, and halves of the four forward "wheels" , can be seen, as can the dense 
array of readout and services cables and the cooling system. The complete MVD 
was installed in ZEUS in April 2001 and has been fully integrated with the ZEUS 
DAQ system; both cosmic-ray and beam-related data have been taken. 

The physics programme addressed by the MVD is that of the flavour decom- 
position of the proton and photon and the search for physics beyond the Standard 
Model. The large increase in luminosity of HERA II, together with the ability to 
tag heavy-quark decays in the MVD, should greatly improve the measurement of 
F% discussed in Section 4.4. After about 500 pb _1 , an uncertainty of around the 



http : / /link . springer .de/link/ service/ journals/10105/ index .html 



EPJdirect Al, 1-11 (2008) Springer- Verlag 



48 



N 1 
> 

0.8 



0.6 



0.4 



0.2 





0.2 0.4 0.6 0.8 1 

K Y 

Fig. 39. Limits on flavour-changing neutral current coupling strength for single 
top production. Limits from HI and ZEUS are plotted on the two-dimensional 
space of the photon and Z coupling strengths. Also shown are similar limits from 
the combined LEP experiments (to the right of the shaded curve) and from CDF 
(to the right of the black shaded lines) . 

2% currently measured on F2 should be obtained. In addition, 6-quark produc- 
tion can be measured precisely; a Monte Carlo simulation J70| of a measurement 
of F 2 b /F 2 c after 500 pb" 1 is shown in Fig. gl]. 

It should also be possible, from a combination of neutral and charged current 
measurements, to separate out the u, d, s, c, b and g contribution to Ft- 

9.1.2 Charged-particle tracking in the forward direction The higher 
luminosity expected at HERA II will increase the number of very high-Q 2 events 
in which the electron or positron is scattered into the forward direction. It will 
also give access to rare processes, including possible physics beyond the Standard 
Model, which tend to have forward jets and/or leptons. The pattern-recognition 
capabilities of the ZEUS Forward Tracker have therefore been improved by the 
replacement of two layers of transition-radiation detector by layers of straw 
tubes. The straws are approximately 7.5 mm in diameter and range in length 
from around 20 cm to just over 1 m. They are constructed from two layers of 50 
/im kapton foil coated with a 0.2 /im layer of aluminium, surrounding a 50 /im 
wire at the centre. The straws are arranged in wedges consisting of three lay- 
ers rotated with respect to each other to give three-dimensional reconstruction. 
Each of the two "supermodules" consists of four layers of such wedges. 
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Fig. 40. A photograph of one half of the MVD, showing the barrel ladders, one 
half of each of the four forward wheels and the cables and services. 

9.1.3 Luminosity monitor The measurement of luminosity at HERA II 
must cope with the greatly increased synchrotron-radiation background and the 
higher probability for multiple bremsstrahlung photons in one beam crossing. To 
compensate for this, two devices, with very different systematic uncertainties, 
have been constructed. Both devices use the information from a small calorime- 
ter placed around 6 m from the interaction point which detects the radiating 
electron. It is hoped that the reduction of systematic error that can be obtained 
from independent luminosity measurements using very different techniques will 
allow a precision of around 1% to be attained. 

9.2 Polarisation 

Polarisations of around 65% have been achieved at HERA I. It is hoped to 
increase the accuracy with which the polarisation can be measured to SP/ P ~ 
2% per bunch per minute. This will be achieved by a collaboration between 
HI, HERMES, ZEUS and the HERA machine in the POL2000 project. The 
collaboration has constructed two instruments, one to measure the longitudinal 
polarisation and the other to measure the transverse polarisation. Both detect 
asymmetries in back-scattered light from high-intensity polarised lasers. 
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Fig. 41. The MC prediction for the ratio of the contribution to F2 of &-quark 
to c-quark production in Q 2 bins as a function of x after 500 pb _1 of data at 
HERA II. 



The combination of high-precision measurements of both luminosity and po- 
larisation will be important in a wide range of HERA II physics, particularly in 
the electroweak sector. The charged current cross section should vanish for the 
appropriate combinations of lepton charge and polarisation. A measurement at 
three polarisations, such as shown in Fig. even with an integrated luminosity 
of only 50 pb _1 per point, will provide an accurate test of this prediction and 
thereby give sensitivity to possible new currents outside the Standard Model. 

Strong polarisation effects are also predicted at high Q 2 in neutral current 
interactions, where, e.g. at Q 2 = 10 4 GeV 2 and x = 0.2, there is a factor of 
two difference between the predicted cross sections for left- and right-handed 
electrons. 

In addition to the use of precise luminosity and polarisation information in 
the study of electroweak processes, polarisation also offers an invaluable tool in 
the study of possible signals beyond the Standard Model. Varying the polarisa- 
tion to reduce the cross sections of Standard Model processes can improve the 
signal to background for new physics signals, such as leptoquarks or supersym- 
metric particles that violate R parity, for which HERA will be competitive with 
the Tevatron for the next few years. 



10 Summary 

We have seen that the study of deep inelastic scattering has been seminal in 
the development of our current understanding of particle physics. In the last 
decade, HERA I has changed our perception of QCD out of all recognition. 
In many cases the precision of the data mandate NNL, or even high order, 
QCD predictions. The study of diffraction and the transition region between 
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Fig. 42. The cross section for charged current interactions. The points at P = 
are obtained from ZEUS preliminary results at the indicated centre-of-mass 
energies, while those at non-zero polarisation are Monte Carlo simulations of the 
expected accuracy in ZEUS assuming the Standard Model cross section for an 
integrated luminosity of 50 pb -1 per point. 



soft and hard physics may be the beginning of the era of quantitative study 
of the central problem of the strong interaction, confinement. The significant 
increase in expected luminosity and the use of new tools such as polarisation 
promise that HERA II physics should continue to provide both new results and 
surprises in the years to come. 
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